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Abstract 

In this paper, we first give a new definition of almost periodic time scales, two 
new definitions of almost periodic functions on time scales and investigate some basic 
properties of them. Then, as an application, by using the fixed point theorem in Banach 
space and the time scale calculus theory, we obtain some sufficient conditions for the 
existence and exponential stability of positive almost periodic solutions for a class of 
Nicholson’s blowflies models on time scales. Finally, we present an illustrative example to 
show the effectiveness of obtained results. Our results show that under a simple condition 
the continuous-time Nicholson’s blowflies models and their discrete-time analogue have 
the same dynamical behaviors. 
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1 Introduction 

To describe the population of the Australian sheep-blowfly and to agree with the ex¬ 
perimental data obtained in |1], Gurney et ah |2] proposed the following delay differential 
equation model: 

x'{t) = —5x{t) + px{t — (1.1) 

where p is the maximum per capita daily egg production rate, 1/a is the size at which the 
blowfly population reproduces at its maximum rate, 6 is the per capita daily adult death rate, 
and r is the generation time. Since equation (11.11) explains Nicholson’s data of blowfly more 
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accurately, the model and its modifications have been now refereed to as Nicholson’s Blowflies 
model. The theory of the Nicholsons blowflies equation has made a remarkable progress in the 
past forty years with main results scattered in numerous research papers. Many important 
results on the qualitative properties of the model such as existence of positive solutions, 
positive periodic or positive almost periodic solutions, persistence, permanence, oscillation 
and stability for the classical Nicholsons model and its generalizations have been established 
in the literature |3l HI El El 13 El 13 113 IHl 1I2]- example, to describe the models of 
marine protected areas and B-cell chronic lymphocytic leukemia dynamics that are examples 
of Nicholson-type delay differential systems, Berezansky et al. [13] and Wang et al. [TT] 
studied the following Nicholson-type delay system: 

{ m 

N[it) = -a,it)N,{t) + /3i(t)iV2(t) + E 

m 

j-i 

where /Sj, c^, 7*^, G C(M, (0, +00)), i = 1, 2, j = 1, 2 ,..., m; in [15], the authors discussed 
some aspects of the global dynamics for a Nicholson’s blowflies model with patch structure 
given by 


x[{t) = -diXiit) + ^ ttijXjit) + ^ /3ijXi{t - Tij)e = 1, 2 ,...,n. 

i=i i=i 

In the real world phenomena, since the almost periodic variation of the environment plays 
a crucial role in many biological and ecological dynamical systems and is more frequent 
and general than the periodic variation of the environment. Hence, the effects of almost 
periodic environment on evolutionary theory have been the object of intensive analysis by 
numerous authors and some of these results for Nicholsons blowflies models can be found in 

[I31I711I31I312312I1I23123]. 

Besides, although most models are described by differential equations, the discrete-time 
models governed by difference equations are more appropriate than the continuous ones when 
the size of the population is rarely small, or the population has non-overlapping generations. 
Hence, it is also important to study the dynamics of discrete-time Nicholson’s blowflies mod¬ 
els. Recently, authors of [231 [25] studied the existence and exponential convergence of almost 
periodic solutions for discrete Nicholson’s blowflies models, respectively. In fact, it is trouble¬ 
some to study the dynamics for discrete and continuous systems respectively, therefore, it is 
significant to study that on time scales, which was initiated by Stefan Hilger (see [26]) in order 
to unify continuous and discrete cases. However, to the best of our knowledge, very few results 
are available on the existence and stability of positive almost periodic solutions for Nichol¬ 
son’s blowflies models on time scales except [27]. But [27] only considered the asymptotical 
stability of the model and the exponential stability is stronger than asymptotical stability 
among different stabilities. 

On the other hand, in order to study the almost periodic dynamic equations on time scales, 
a concept of almost periodic time scales was proposed in [2S]- Based on this concept, almost 
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periodic functions [28], pseudo almost periodic functions [29], almost automorphic functions 
[30] , weighted pseudo almost automorphic functions [31], weighted piecewise pseudo almost 
automorphic functions [32] and almost periodic set-valued functions [33] on on time scales 
were dehned successively. Also, some works have been done under the concept of almost 
periodic time scales (see [311 1351 136] 133 EH] 139] 001 SI])- Although the concept of almost 
periodic time scales in [28] can unify the continuous and discrete situations effectively, it is 
very restrictive. This excludes many interesting time scales. Therefore, it is a challenging and 
important problem in theories and applications to hnd new concepts of periodic time scales 

([^ 031 03 03 iS]). 

Motivated by the above discussion, our main purpose of this paper is hrstly to propose a 
new dehnition of almost periodic time scales, two new dehnitions of almost periodic functions 
on time scales and study some basic properties of them. Then, as an application, we study 
the existence and global exponential stability of positive almost periodic solutions for the 
following Nicholson’s blowflies model with patch structure and multiple time-varying delays 
on time scales: 

n 

xf{t) = -Ci{t)xi{t) + ^ bik{t)xk{t) 

k=l,k^i 

n 

+ ^(^ij{t)xi{t - i = 1, 2,..., n, (1.2) 

i=i 

where f G T, T is an almost periodic time scale, Xi{t) denotes the density of the species in 
patch i, bik{k ^ i) is the migration coefficient from patch k to patch i and the natural growth 
in each patch is of Nicholson-type. 

For convenience, for a positive almost periodic function / : T —)■ M, we denote /■*■ = 
sup/(t), f~ = inf fit). Due to the biological meaning of fll.2p . we just consider the following 

teT 

initial condition: 

(Pi{s) >0, s e [to - 6, to]T, to e T, i = 1, 2,..., n, (1.3) 

where 9 = maxsup{rij(t)}, [to — 9, toji = [to — 9, to] H T. 
h-i) teT 

This paper is organized as follows: In Section 2, we introduce some notations and deh¬ 
nitions which are needed in later sections. In Section 3, we give a new dehnition of almost 
periodic time scales and two new dehnitions of almost periodic functions on time scales, and 
we state and prove some basic properties of them. In Section 4, we establish some sufficient 
conditions for the existence and exponential stability of positive almost periodic solutions of 
fll.21) . In Section 5, we give an example to illustrate the feasibility of our results obtained in 
previous sections. We draw a conclusion in Section 6. 

2 Preliminaries 

In this section, we shall hrst recall some dehnitions and state some results which are used 
in what follows. 
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Let T be a nonempty closed subset (time scale) of M. The forward and backward jump 
operators cr, p : T —)■ T and the graininess p : T —)■ M+ are dehned, respectively, by 

a{t) = inf{s G T : s > f}, p{t) = sup{s G T : s < t} and p(t) = a{t) — t. 

A point t G T is called left-dense if t > infT and p{t) = t, left-scattered if p{t) < t, 
right-dense if t < supT and a{t) = t, and right-scattered if a{t) > t. If T has a left-scattered 
maximum m, then = T \ {m}; otherwise = T. If T has a right-scattered minimum m, 
then Tfc = T \ {m}; otherwise = T. 

A function / : T —)■ R is right-dense continuous provided it is continuous at right-dense 
point in T and its left-side limits exist at left-dense points in T. If / is continuous at each 
right-dense point and each left-dense point, then / is said to be continuous function on T. 

For p : T —)• R and t G we dehne the delta derivative of y{t), to be the number 

(if it exists) with the property that for a given £ > 0 , there exists a neighborhood 17 of t such 
that 

\[yia{t)) - 2 /(s)] - y^{t)[a{t) - s]| < e\a{t) - s\ 

for all s E U. 

If y is continuous, then y is right-dense continuous, and if y is delta differentiable at t, 
then y is continuous at t. 

Let y be right-dense continuous. If Y^(t) = y{t), then we dehne the delta integral by 
vis) As = Y{t)-Y{a). 

A function r : T —)■ R is called regressive if 1 -|- p{t)r{t) ^ 0 for all t G T^. The set of all 
regressive and rd-continuous functions r ; T —)■ R will be denoted hj 71 = 7?.(T) = 7^(T, R). 
We dehne the set Tl'^ = 7^+(T, R) = {r G 7?. ; 1 -|- p{t)r{t) >0, Vf G T}. 

Lemma 2.1. ([^) Suppose that p G , then 
(i) ep{t, s) > 0, for all t,s E T; 

{ii) if pit) < q(t) for all t > s,t, s E T, then Cpit, s) < eq(t, s) for all t > s. 

Definition 2.1. A subset S o/R is called relatively dense if there exists a positive number 
L such that [a, a -|- L] fi S' 7 ^ 0 for all a G R. The number L is called the inclusion length. 

Definition 2.2. A time scale T is called an almost periodic time scale if 

n={rGR:t±rG T,Vf G T} ^ {0}. 

The following dehnition is a slightly modihed version of Dehnition 3.10 in [28] . 

Definition 2.3. Let T be an almost periodic time scale. A function f E ^(T x i7,E”) is 
called an almost periodic function in t eT uniformly for x E D if the e-translation set of f 

E{e,f,S} = {r G n : \ f{t + r,x) - f{t,x)\ < e, V(t,x) G T x S'} 

is relatively dense for all e > 0 and for each compact subset S of D; that is, for any given 
£ > 0 and each compact subset S of D, there exists a constant l{e,S) > 0 such that each 
interval of length S') contains a T{e, S) G E{e, f,S} such that 

\f{t + T,x) - fit, x)\ < e, V(f,x) G T X S'. 

r is called the e-translation number of f. 
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3 Almost periodic time scales and almost periodic func¬ 
tions on time scales 

In this section, we will give a new definition of almost periodic time scales and two new 
definitions of almost periodic functions on time scales, and we will investigate some basic 
properties of them. Our new definition of almost periodic time scales is as follows: 

Definition 3.1. A time scale T is called an almost periodic time scale if the set 

n:={reT:T, ^0 and T, ^ {0}} ^ {0}, 

where T.r = Tn{T — r} = Tn{f — r:fG T}, satisfies 

(t) n ^ 0, 

(a) if Ti, T 2 G n, then Ti ± r 2 G 11, 

{Hi) f := T(n) = n T, ^ 0. 

rsn 

Clearly, if f G then t + r G T. If f G T, then f + r G T for r G 11. 

Remark 3.1. Obviously, if T is an almost periodic time scale under Definition \3.1[ then 
infT = —oo and supT = +cxo. If T is an almost periodic time scale under Definition 12.21 
then T is also an almost periodic time scale under Definition \S.l\ and in this case, T = T. 

Example 3.1. Let T = Z U {|}. For every r ^ X, we have T,- = Z and Ti = {0}. Hence 

n = Z and T = Hrsn^ 7 ^ 0- So, T is an almost periodic time scale under Definition 
\3.1\ but it is not an almost periodic time scale under Definition\2fM 

Lemma 3.1. IfT is an almost periodic time scales under Definition AS.li then T is an almost 
periodic time scale under Definition \2.S\ 

Proof. By contradiction, suppose that there exists a. to E T such that for every r G 11 \ {0}, 

to + T or to — T ^T. ^ 

Case (i) li to + t ^ T, then there exists a rtp G 11 such that to + t ^ On one 

hand, since to + t E T, to + t + ^ T. On the other hand, since to ^ T and t + E 11, 

to + T + Tto G T. This is a contradiction. 

Case (ii) to — r then there exists a ft^ E Ti such that to — r ^ On one 

hand, since to — t E T , to — r + ^ T. On the other hand, since H E T and — r + ft^ G II, 
to — T + ftQ G T. This is a contradiction. 

Therefore, for every f G T, there exists a r G II \ {0} such that f ± r G T. Hence, T is an 
almost periodic time scale under Definition 12.21 The proof is complete. | 
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Throughout this section, E"" denotes M” or C"", D denotes an open set in E” or D = E", 
S denotes an arbitrary compact subset of D. 

From [2B], under Dehnitions 12.21 and 12.31 we know that if we denote by BUC{T x D,W^) 
the collection of all bounded uniformly continuous functions from T x D to then 

AP{T X D, M") C BUC{T x D, M"), (3.1) 

where AP{T x D,W^) are the collection of all almost periodic functions in t G T uniformly 
for X ^ D. It is well known that if we let T = M or Z, fl3.1|) is valid. So, for simplicity, we 
give the following dehnition: 

Definition 3.2. Let T be an almost periodic time scale under sense of Definition \d.l[ A 
function f G BUC{T x D,E”') is called an almost periodic function in f G T uniformly for 
X E D if the e-translation set of f 

E{e,f,S} = {r G n ; \ f{t + T,x) - f{t,x)\ < e, W{t,x) G T x 5} 

is relatively dense for all e > 0 and for each compact subset S of D; that is, for any given 
£ > 0 and each compact subset S of D, there exists a constant l{e,S) > 0 such that each 
interval of length S) contains a T{e, S) G E{e, /, S} such that 

\f{t + T,x)- f{t,x)\<e, \/{t,x)ETxS. 

This T is called the e-translation number of f. 

Remark 3.2. If T = M, then T = R, in this case, Definition 13.31 is actually eguivalent to 
the definition of the uniformly almost periodic functions in Ref. // T = Z, then T = Z, 

in this case. Definition 13.SI is actually eguivalent to the definition of the uniformly almost 
periodic seguences in Refs. 

For convenience, we denote by AP{T x D,EL) the set of all functions that are almost 
periodic in t uniformly for x E D and denote by AP{T) the set of all functions that are 
almost periodic in t G T, and introduce some notations: Let a = {an} and (3 = {(3n} be two 
sequences. Then (3 C a means that /? is a subsequence of a; a -\- (3 = + (3n}] —o. = {—an}; 

and a and (3 are common subsequences of a and (3 , respectively, means that 
and [3n = (3'n(^-^ for some given function n{k). We introduce the translation operator T, 
Tafit, x) = git, x) means that g{t, x) = lim f{t + x) and is written only when the limit 

n—^-\-oo 

exists. The mode of convergence, e.g. pointwise, uniform, etc., will be specihed at each use 
of the symbol. 

Similar to the proofs of Theorem 3.14, Theorem 3.21 and Theorem 3.22 in [22], respectively, 
one can prove the following three theorems. 

Theorem 3.1. Let f E UBC{T x D,¥L), if for any seguence a C 11, there exists a (Z a 
such that Taf exists uniformly onT x S, then f E AP{T x D,'E'^). 
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Theorem 3.2. If f & AP{T x then for any e > t), there exists a positive constant 

L = L{e, S), for any a G M, there exist a constant rj > 0 and a G M such that ([a, a+rj] fill) C 
[a, a + L] and ([a, a + r;] n II) C E{e, /, S). 

Theorem 3.3. If f, g G AP{TxD, E"-), then for any £ > 0, E{f, e, S)nE{g, e, S) is nonempty 
relatively dense. 

According to Definition 13.21 one can easily prove 

Theorem 3.4. If f & AP{T x D, E""), then for any a G M, 6 G 11, functions af, f{t + b,-) G 
AP{T X D,E^). 

Similar to the proofs of Theorem 3.24, Theorem 3.27, Theorem 3.28 and Theorem 3.29 in 
[25] . respectively, one can prove the following four theorems. 

Theorem 3.5. If f,g G AP(T x D, E"), then f + g, fg ^ AP(T x D, E""), if inf \g(t, x) | >0, 
then f/g G AP(T x D,E"'). 

Theorem 3.6. If fn G AP{T X D,W^){n = 1,2,...) and the sequence {fn} uniformly con¬ 
verges to / on T X S', then f G AP(T x D, E"). 

Theorem 3.7. Iff G AP(TxD,E^), denote Eft,x) = f(s, x)As, then E G AP(TxD,K^) 
if and only if E is bounded on T x S. 

Theorem 3.8. If f & APfl X D,E"'), E{-) is uniformly continuous on the value field of f, 
then E o f is almost periodic in t uniformly for x & D. 

By Definition 13.21 one can easily prove 

Theorem 3.9. Let / : M —)■ M satisfies Lipschitz condition and (p{t) G AP{T), then f{ip{t)) G 
AP{T). 

Definition 3.3. Let Aft) be an n x n rd-continuous matrix on T, the linear system 

x‘^{t) = A{t)x{t), f G T (3.2) 

is said to admit an exponential dichotomy on T if there exist positive constant k, a, projection 
P, and the fundamental solution matrix X ft) of fl3.2p . satisfying 

\X{t)PX-^{a{s))\ < keQ^{t,a{s)), s,t e T, t > a(s), 


jX(t)(I - P)X \a(s))j < kee^(o-(s),t), s,tET, t < a(s), 
where | • | is a matrix norm on T, that is, if A = {aij)nxm, then we can take |A| = 

n m 1 

(EE kit)’. 

1=1 j=l 

Similar to the proof of Lemma 2.15 in |12], one can easily show that 
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Lemma 3.2. Let aii{t) be an uniformly bounded rd-continuous function on T, where aii{t) > 
0, —Quit) G for every t G T and 

min {inf aait)} > 0, 

l<i<n t&T 

then the linear system 

x^{t) = diag{-aii{t), -a 22 {t),-ann{t))x{t) 
admits an exponential dichotomy on T. 

According to Lemma 13.li T is an almost periodic time scales under Definition 12.21 we 
denote the forward and the backward jump operators of T by a and p, respectively. 

Lemma 3.3. If t is a right-dense point on T, then t is also a right-dense point on T. 

Proof. Let f be a right-dense point on T, then 

t = a{t) = inf{s G T : s > f} > inf{s G T : s > f} = cr(f). 

Since a(t) >t,t = o'(t). The proof is complete. | 

Similar to the proof of Lemma 13.31 one can prove the following lemma. 

Lemma 3.4. If t is a left-dense point T, then t is also a left-dense point on T. 

For each / G C(T,R), we define / : T —)■ R by f{t) = f{t) for f G T. From Lemmas 13.31 
and 13.41 we can get that / G ^(T, R). Therefore, F dehned by 

F{t) := [ fir)A t, to,t ef 
dto 

is an antiderivative of / on T, where A denotes the A-derivative on T. 

Set n = {rGn;f±rG T}. We give our second definition of almost periodic functions 
on time scales as follows. 

Definition 3.4. Let T be an almost periodic time scale under sense of Definition \3.1[ A 
function f G BUC{T x D,K^) is called an almost periodic function in t G T uniformly for 
X E D if the e-translation set of f 

E{e,f,S} = {r G n ; \ f{t + T,x) - f{t,x)\ < e, V(f,x) G T x S'} 

is relatively dense for all e > 0 and for each compact subset S of D; that is, for any given 
e > 0 and each compact subset S of D, there exists a constant l{e,S) > 0 such that each 
interval of length l{e, S) contains a T{e, S) G E{e, f, S'} such that 

\f{t + T,x)- f{t,x)\<e, V(f,x)GfxS'. 


This T is called the e-translation number of f. 









Remark 3.3. It is clear that if a function is an almost periodic function under Definition 
\3.^A then it is also an almost periodic function under Definition \3.Ji\ 

Remark 3.4. SinceT is an almost periodie time scales under Definition \2.3[ under Definition 
\3.^ all the results obtained in fMBf remain valid when we restrict our discussion to T. 

In the following, we restrict our discuss under Definition 13.41 
Consider the following almost periodic system: 

x^{t) = A{t)x{t) + fit), t G T, (3.3) 

where Ait) is a n x n almost periodic matrix function, fit) is a n-dimensional almost periodic 
vector function. 

Similar to Lemma 2.13 in |12], one can easily get 

Lemma 3.5. If linear system fl3.2p admits an exponential dichotomy, then system fl3.3p has 
a bounded solution xit) as follows: 

/ t /* + c>o 

Xit)PX-\ais))fis)As - / Xit)il - P)X-\ais))fis)As, t e T, 

-oo J t 

where Xit) is the fundamental solution matrix of (13.2p . 

By Theorem 4.19 in [28], we have 

Lemma 3.6. Let Ait) be an almost periodic matrix function and fit) be an almost periodic 
vector function. If fl3.2l) admits an exponential dichotomy, then fl3.3|) has a unique almost 
periodic solution: 

/ i ^ ^ ^ ^ r+co ^ ^ ^ ^ ^ 

Xit)PX-\ais))fis)As - / Xit)il - P)X-\ais))fis)As, t e T, 

-OO J t 

where Xit) is the restriction of the fundamental solution matrix of fl3.2p on T. 

From Dehnition 13.21 and Lemmas 13.51 and 13.61 one can easily get the following lemma. 

Lemma 3.7. If linear system fl3.2p admits an exponential dichotomy, then system 03.31) has 
an almost periodic solution xit) can be expressed as: 

/ t r-\-oo 

Xit)PX-\ais))fis)As - / Xit)il - P)X-\ais))fis)As, t e T, 

OO J t 

where Xit) is the fundamental solution matrix of 03.2p . 
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4 Positive almost periodic solutions for Nicholson’s blowflies 
models 


In this section, we will state and prove the sufficient conditions for the existence and 
exponential stability of positive almost periodic solutions of fll.2p . Throughout this section, 
we restrict our discussion under Dehnition 13.41 

Set B = {93 G C(T, M”) : ip = [ipi, ip<ii ..., ipn) is an almost periodic function on T} with the 
norm ||</3 ||b = max sup \ipi{t)\, then B is a Banach space. Denote C = C{[to — 0 ,^ 0 ]?,®”) and 

l<i<n jgT 

C{Ai, A 2 } = {(p = {(pi, (P 2 , ■ ■ ■ ,^n) ^ c : Ai < ipi{s) < A 2 ,s e [to - 0, to]T, i = 1,2, ... ,n}, 
where 0 < Ai < A 2 are constants. 

In the proofs of our results of this section, we need the following facts: There exists a 
unique G (0,1) such that ^ ~ 0.7215354) and sup The function xe~^ 

x><; 

decreases on [l,+oo). 

Lemma 4.1. Assume that the following conditions hold. 

(Hi) Ci,bik,/3ij,aij,Tij G AP(T,M+) and c~ > > 0, t - Tij{t) G T, 

i,k,j = l,2,...,n. 

{H 2 ) t §<lD = l,2,...,n. 

(Ps) There exist positive constants ^1,^42 satisfy 


A 2 > max 

l<2<n 


1 



k=l^k^i ^ - 



and 


mm 

l<2<n 


k=l,k^i ^ j=l ^ 


> Ai> 


min 




Then the solution x{t) = [x\{t), X 2 {f),... ,Xn{f)) of (11.21) with the initial value ip G C{Ai,A 2 } 
satisfies 


< Xi{t) < ^ 2 , t e [to, +oo)t, i = 1, 2,..., n. 

Proof. Let x{t) = x{t;to,ip), where ip G C{y4i,y42}. At hrst, we prove that 

Xi{t) < A 2 , t G [to, ^{ip))^, i = l,2,...,n, (4.1) 

where [to,'r]{ip))j is the maximal right-interval of existence of x{t-,to,ip). To prove this claim, 
we show that for any p > 1, the following inequality holds 

Xi{t) < pA 2 , t G [to, h(<p))T, i = 1,2,... ,n. (4.2) 
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By way of contradiction, assume that fl4.2p does not hold. Then, there exists io ^ {1, 2,..., n) 
and the hrst time ti G such that 

Xioiti) > pA2, Xi^(t) < pA2, t e [to - 

Xk{t) < pA 2 , for t e [to - 9, tijx, k = 1,2,... ,n. 

Therefore, there must be a positive constant a > 1 such that 

Xioih) = apA2, Xi^{t) < apA2, t e [to - 9,ti)T, 

Xk{t) < apA 2 , ioik ^ i^, t G [to - 6',ti]T, A; = 1, 2,..., n. 

In view of the fact that supMe““ = - and ap > 1, we can obtain 

n >0 ^ 


= -Cio(ti)xio(ti) + ^ bi^k{ti)xk{ti) 

k=l,k^iQ 




< -Ci,apA2+ Y1 Kk(^P^2 + Yl 

k=l,k^iQ 


A 


loj 


j=l ^ioj ^ 


< apa 


*0 


A 2 + ^ 


n 4 L+ n 






C- ' ^ C CY- -6 

k=l,k^iQ ^0 j=l ^0 


< 0 , 


which is a contradiction and hence fl4.2p holds. Let p —)■ 1, we have that fl4.ip is true. Next, 
we show that 


Xi{t) > Ai, t G [to, p(<p))t, t = 1, 2, ..., n. (4.3) 

To prove this claim, we show that for any / < 1, the following inequality holds 

Xi{t) >lAi, t G [to, p((p))t, i = 1, 2,..., n. (4.4) 

By way of contradiction, assume that fl4.4p does not hold. Then, there exists A G {1, 2,... , n} 
and the hrst time t 2 G [to,p{ip))j such that 

Xi^(t 2 ) < lAi, Xi^(t) > I, t G [to - 0,t2)T, 

Xk{t) > lAi, for k ii, t G [to - 9, t 2 ]T, A: = 1, 2,..., n. 

Therefore, there must be a positive constant c < 1 such that 

Xi^{t 2 ) = clAi, Xi^{t) > cl, t G [to - 6 ',t 2 )T, 

Xkit) > clAi, for k ^ ii, t G [to — 9, t 2 ]T, A; = 1, 2,..., n. 
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Noticing that c/ < 1, it follows that 


0 > Kk{h)Xk{t2) 

k=l,k^ii 

n 

i=i 


> -c+cMi + ^ b-^i^clAi + A 




i=i 

n 


a 


*1J 


dc;(-/ii+ 5^ + 

/C = l,fc^ 2 l 


which is a contradiction and hence (14.41) holds. Let / —>■ 1, we have that (14. 3 h is true. Similar 
to the proof of Theorem 2.3.1 in [51j, we easily obtain ri{^p) = +cx3. This completes the proof. 


Remark 4.1. IfT = M, then jji{t) = 0, so, —Ci G TA. //T = Z, then jn{t) = 1, so, —Ci G TA 
if and only if Ci < 1. 

Theorem 4.1. Assume that {Hi) and (R 3 ) hold. Suppose further that 

(R 4 ) —Ci G TZ^, where VA denotes the set of positive regressive functions, i = 1,2,... ,n. 

n ^11+ 

(^ 5 ) E Kk + J2^<cfA = ^A,---,ri. 

k=l^k^i j=l 

Then system (11.2P has a positive almost periodic solution in the region B* = |(p| (p G B, Ai < 
Proof. For any given <^9 G B, we consider the following almost periodic dynamic system: 

n 

x^{t) = -Ci{t)xi{t) + ^ bik{t)(pk{t) 

k=l.,k^i 

n 

+ /3At)Ti{t - r,,(t))e-“-W), t = l,2,...,n. (4.5) 

i=i 

Since mini<j<„{c)"} > 0, f G T, it follows from Lemma [3.21 that the linear system 

xf (f) = -Ci{t)xi{t), i = 1,2,... ,n 

admits an exponential dichotomy on T. Thus, by Lemma 13.71 we obtain that system fl4.5p 
has an almost periodic solution x^ = {x^.^,x^p^, ..., x^^), where 
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As, i = 1, 2,..., n. 


j=i 

Define a mapping T : B* —)■ B* by 

T(f{t) = x^pit), 'iip e B*. 

Obviously, B* = {ip\ 99 G B, Ai < (pi{t) < A 2 ,t G T, i = 1, 2,..., n} is a closed subset of B. 
For any (p G B*, by use of {H 2 ), we have 


XpAt) < / e_^-{t,a{s)) 


n n 0+ ■ 

^ k=l,kpi j=l 


As 


< 


C: 


Y1 + 

k=l,kpi 


Atj 1 

^ X - 

■ 1 “i7 ^ 

J=l V 


< A 2 , f = 1, 2,... ,n 


and we also have 


Xp.{t) > / e_^+(f,a(s)) 


^k=l,kpi j=l 


As 


> 


^ A 2 / 3 , 

k=l,kpi j=l 




> Ai, i = 1, 2,..., n. 

Therefore, the mapping T is a self-mapping from B* to B*. 

Next, we prove that the mapping T is a contraction mapping on B*. Since sup | 


l-u l _ J_ 

pu \ p 2 ? 


we hud that 


u><; 


\xe ^ — ye '^\ = 


1 - (a^ + ^(l/-a:)) 


^x+^{y-x) 


lx - y\ 


< -^\x-y\, x,y ><;, 0 < ^ < 1. 

For any p = {pi, p2 ,..., Pn)'^, A = (t/’n t/’2, • • •, An)'^ £ we obtain that 
|(Tv.),(t)-(T^),(t)| 


< 


bik{s){pk{s) - ^/’fc(s))As 

k=l,kpi 

r a(s)) 0.(s " r,j(s))e-“«W«’-<*-’-'''‘» 


i=i 


-^i(s - Ti,{s))e-^^As)i’i{s-nps))\ 


As 
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1 ^ ^ iQ / \ 

< — 5^ V’||b+ / e_c,(t,cT(s))^^^(aij(s)(^i(s-ry(s)) 

i=i ^ 

/I ” ” /3+ \ 

< l;:^ ^+fc + 5^;;^)||v?-^/^llB, * = l,2,...,n. 


* k=l,k^i 


■ 1 Cj e2 

j=i * 


It follows that 

||T0-T^||b< II^-V'IIb, 

which implies that T is a contraction. By the hxed point theorem in Banach space, T has 
a nniqne hxed point ip* G B* snch that Tip* = ip*. In view of fl4.5p . we see that ip* is a 
solntion of fll.2l) . Therefore, fll.2p has a positive almost periodic solution in the region B*. 
This completes the proof. | 


Definition 4.1. Let x*{t) = {xl{t), xl{t),... ,x*^{t)Y be an almost periodic solution of fll.2p 
with initial value ip*{s) = {ipl{s), ip^^s),..., ip'^{s))'^ G C'{Ai,A 2 }. If there exist positive 
constants A with ©A G TZ^ and M > 1 such that such that for an arbitrary solution x(t) = 
{xi{t),X 2 {t),... ,Xn{f)Y of 01.2|) with initial value ip{s) = {ipi{s), ip 2 {s ),..., ipn{s)Y' G C{Ai, A 2 } 
satisfies 

llx - x*|| < M\\ip - ip*\\eQxit,to), to G [-9, oo)t, t > to, 
where ||(p —(p*||o = max < sup \‘Pi{t) ~ Tl{t)\ [■ for ip,ip G C{Ai, A 2 }■ Then the solution 

i<i<n [ te[to-e,to] J 

x*(t) is said to be exponentially stable. 

Theorem 4.2. Assume that {Hi), {H 3 )-{H 5 ) hold. Then the positive almost periodic solution 
x*{t) in the region B* of 01.2p is unique and exponentially stable. 

Proof. By Theorem 14.11 01.2p has a positive almost periodic solution x*{t) in the region 
B*. Let x{t) = {xi{t), X 2 {t),... ,Xn{t))^ be any arbitrary solution of 01.21) with initial value 
ip{s) = ((pi(s), ip 2 {s), • • •, ipn{s))'^ G C{Ai, A 2 }. Then it follows from 01.2p that for t >to,i = 
1,2,... ,n. 


{x,{t) - x*{t))^ 

n 

= -Ci {t) {Xi {t) - X* (f)) + ^ bik (t) {Xk {t)-x*k (t)) 

n 

+ Pijit) [xi{t - - x*{t - . (4.6) 

i=i 

The initial condition of 04.6p is 

A{s) = P>i{s) - x*{s), s G [to - d,to]T, i = l,2,...,n. 
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For convenience, we denote Ui(t) = Xi(t) — x*(t),i = 1,2,... ,n. Then, by fl4.6p . we have 


Ui{t) = Wi(to)e_ci(t,to) + / e_c,(t,a(s)) ^ bik{s)uk{s)As 

“^*0 k=l,k^i 


+ / e_c,(t,cr(s))^/3y( 
Jto 


-x*{s-Tij{s))e 


-aij{s)x*(s-Tij(s)) 


As, t > to,i = 1,2,... ,n. 


(4.7) 


For o; G M, let Fj(ci;) be defined by 


sST 


Fi(a;) = Ci -a;-exp{a;sup/x(s)}( ^ 6,:^ + — ^/3j expjwrj} ), i = 1, 2,..., n. 


k=l,k^i 


i=i 


In view of {H 2 ), we have that 


r 4 ( 0 ) = c, -( ^ 6++) > 0,i = l,2,...,n. 


k=l,k^i 


J=1 


Since Fj(a;) is continuous on [0, + 00 ) and Fj(a;) —00 as u ^ + 00 , so there exists cuj > 0 
such that Fj(ci;j) = 0 and Fi(a;) > 0 for a; G {0,oji),i = 1,2,... ,n. By choosing a positive 
constant a = min \^ui,uj 2 ,... ,ujn], we have Fj(a) > 0, i = 1, 2,..., n. Hence, we can choose a 
positive constant 0 < a < min |a, min {c“}} such that 


\<i<n 

Fi(a) >0, i = 1,2,... ,n. 


which implies that 

exp{Q; sup /i(s)} 


sGT 


C: - a 


^tk + expjarj} ) <1, i = 1, 2,..., n. 


Take 


k=l^k^i 


M = max 

l<2<n 


j=l 


Y1 Kk + h 12, 




k=l^k^i j=l 

It follows from {H^) that M > 1. Besides, we can obtain that 

exp{a sup/i(s)} 


M 


< 


sGT 


C: - a 




■ = jl = l 

In addition, noticing that eQa{t,to) > 1 for f G [to — 9,to]T. Hence, it is obvious that 

||m||b < Mee„(t, to)Iloilo, Vt G [to - 6',to]T- 
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We claim that 


(4.8) 


II'uIIb < Mee„((,(o)||i;i||o, V( e ((o,+oo)t. 

To prove this claim, we show that for any p > 1, the following inequality holds 


mIIb < pMee„(«, io)||^||i,, V(6 (i|,,+oo)T, 


(4.9) 


which implies that, for i = 1,2,...,n, we have 


Ui{t)\ < pMeQa(,t,to)\\'ip\\o, \ft e (fo,+oo)T. 


(4.10) 


By way of contradiction, assume that fld.lOp is not true. Then there exists ti G (to, +oo)t and 
Zq G {1, 2,..., 77,} such that 

\uio{ti)\ > pMeea(fi,fo)||'0||o, \uio{i)\ < to)||'0||o, te (to,^i)T, 

\uk{t)\<pMeQa{t,to)\\^jJ\\o, ioik^io, f G (to, ^i]t, ^ = 1, 2,..., n. 

Therefore, there must be a constant q > 1 such that 

l“*o(^i)l = QpMeea{tl,to)\\^|J\\o, \uio{t)\ < gpMeea(t,to)||'0l|o, t e (to,^i)T, 

\uk{t)\ < qpMeQa{tl,to)\\^p\\o, t E k = 1,2,... ,n. 

According to fl4.7p . we have 









< e_ci^{ti,to)\\'ip\\o + qpMeQo,{ti,to)\\'ip\\o 
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,=i e 

= qpMeQa{tuto)\\i>\\o 


o + gpMee„(ti,to)||V'l|oexp{asup/i(s)}( V 6+^ 

sGT ' ^' 


k=l,k^iQ 


qpM 


e-ciQffia(^i,^o)+exp{asupp(s)} V 6+ 


m /3+ \ /.*] 

£^^exp{ari}) / 

j_l / J to 


e-Ci^®a{tu(T{s))As 


< gpMeea(ti,to)||'0l|o 


A 




j=i 

< gpMeea(ti,i:o)||Al|o< 


ICio - « 


-«)(^iAo) + exp{asupp(s)} E '’.V 

^k=l,k^io 

rti 

) Jt ~ a(s))As 


exp{asupp(s)} 

1 sGT 


qpM 




EA^='p{“A> 


i=i 


e_(c-_«)(^i,^o) + 

^ in > 


E 

fc = l,fc:^ 2 Q 

exp{asup /i(s)} 

sGT 




^ ^ ^iofc 
/c=l,fc:7^io 


^/3Aexp{arE} 


i=i 

< gpMeea(ti,to)||Al|o< 


exp{asup p(s)} 

1 _ sGT 

M 




5^Aiexp{arE} 


i=i 


e_(c-_a)Al,^o) + 

V IQ ^ 


Kk 

k=l,k^io 

exp{a sup p(s)} 

sGT 


Cio -« 


E 

k=l,k^io 


5^/3+.exp{arE} 


i=i 


< gpMeeaAi,to)||Al|o, 


which is a contradiction. Therefore, fld.lUp and fl4.9p hold. Let p ^ 1, then fl4.8p holds. 
Hence, we have that 

||m||b < M\\'il!\\oeQait,to), t e [to,+ oo)t, 

which implies that the positive almost periodic solution x*{t) of fll.2p is exponentially stable. 
The exponential stability of x*{t) implies that the uniqueness of the positive almost periodic 
solution. The proof is complete. | 
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Remark 4.2. It is easy to see that under definitions of almost periodic time scales and almost 
periodic functions in the conclusions of Theorems \f.l\ and \f.^ are true. 


Remark 4.3. From Remark \4.1\ Theorem \4.1\ and Theorem 14-21 we can easily see that 
if Ci{t) < = 1,2,... ,n, then the continuous-time Nicholson’s blowflies models and the 

discrete-time analogue have the same dynamical behaviors. This fact provides a theoretical 
basis for the numerical simulation of continuous-time Nicholson’s blowflies models. 

Remark 4.4. Our results and methods of this paper are different from those in |F7| /. 

Remark 4.5. When T = M or T = Z, our results of this section are also new. If we take 
T = M, y4i = 1, ^42 = e, then Lemma \4.1l Theorem \4.1\ and Theorem \4.1\ improve Lemma 2-4, 
Theorem 3.1 and Theorem 3.2 in respectively. 

5 An example 

In this section, we present an example to illustrate the feasibility of our results obtained 
in previous sections. 

Example 5.1. In system (11.211 . let n = 3 and take coefficients as follows: 

Ci{t) = 0.21 + 0.01 sin ^ buit) = 0.03 + 0.01 cosvrt, bifit) = 0.06 + 0.01 cos Vst, 

= 0.07 + 0.02sinvrt, fii 2 {t) =0.15 + 0.01 cos\/3t, 
fiisit) = 0.15 + 0.01 sin y = 0 ^ 12 ( 4 ) = aifit) = 0.91 + 0.09| sin V3t\, 
rn(t) = Tuit) = e°-^l5)1,ri3(t) = 

C 2 {t) = 0.3 + 0.02 sin ; ^21 (^) = 0.05 + 0.01 cos Vst, h 22 ,{t) = 0.05 + 0.01 sin V2t, 

fi 2 i{t) = 0.06 + 0.01 cosvrt, /322(^) = 0.04 + 0.01 cos \/3f, 1423,(4) = 0.09 + 0.01 cos ^g^^ > 

ct 2 i{t) = 0.8 + 0.2 sin y/2t, a 22 {t) = 0.8 + 0.2 cos V^t, a 23 {t) = 0.8 + 0.2 sinTTf, 
r2i(t) = e°-2|“")"‘+5)l,r22(t) = r23(t) = 

^ 3 ( 4 ) = 0.41 + 0.01 sin ^g^^ ! ^31 (^) = 0.16 + 0.01 sin Vst, 132 ( 4 ) = 0.13 + 0.01 cos \/2t, 

= 0.02 + 0.01 cos , (432{t) = 0.032 + 0.01 cos V2t, 

(I^fit) = 0.022+ 0.001 sin , a 3 i(t) = 0.8 + 0.2| sin \/3t|. 
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C( 32 (t) = 0.8 + 0.2 sin Vst, Q: 33 (t) = 0.8 + 0.2 sin ) 

r3i(t) = r33(t) = 

By calculating, we have 

cY = 0.2 ,4 = 0.22, bY 2 = 0.02, 6 + = 0.04,6^3 = 0.05, = 0.07, 

/3n = 0 . 05 ,/ 3 + = 0.09,4-2 = 0 . 14 ,/ 3 + = 0.16,/3Y3 = 0.14,/3+ = 0 . 16 , 
a^i = a^2 ~ (^13 ~ 0 . 91 , = 0:^2 ~ ~ 

4 = 0.28, 4 = 0.32, bYi = 0.04, b^i = 0.06, 6^3 = 0.04, 6^3 = 0.06, 
/3Y, = 0.05,13+ = 0.07,13^2 = 0.03,/3+ = 0.05,13^3 = 0.08,13+ = 0.1, 

*^21 ~ 0 . 6 , a 2 i = 1 , ^22 ~ 0 . 6 , 0:^2 ~ 1 ) ^23 ~ 0 . 6 , 0-23 = 1 , 

c- = 0.4, C 3 + = 0.43, 63-1 = 0.15, 63 + = 0.17, bY 2 = 12, b+ = 0.14, 

4 -, = 0.01,/3+ = 0 . 03 , 4-2 = 0.022,/3+2 = 0.042,47, = 0.21,/3+3 = 0.23, 

0:31 = 0.8, o 4 i — 1) 0^32 ~ 0.6, a ^2 ~ 1) CI33 ~ 0.6, = 1. 


Hence, 


^ Kk ^ 


12 


13 


_ I _ 


k=l,k7^1 1 


0.04 + 0.07 
02 


= 0.55 < 1, 


+ + ^ ^ 0.06 + 0.06 ^ 
c,- cj cj 


0.28 


^ 3 fc ^31 ^32 0.17 + 0.14 

> -Ih = ^ ^ ^-= 0.775 < 1, 

n- n- n- ^ ’ 


A:=l,fc^3 3 


+ 6 ^. + 4■ + ^ + ^ 


^12 T 1^13 


= 0.04 + 0.07 + 


0.4 


0.09 0.16 0.16 


= 0.1655 < ci 
e“ 

^ . 0.06 + 0,06 + ^ + 1++ 1+= 0.1457 < c; 


1^31 I 1^32 I /^33 


bti + bt. +4++ 47+ = 0.03 + 0.04 + 


^32 


e 2 

+ 2 

+ 2 

0.07 

0.05 

0.07 

e 2 

+ 2 

+ 2 

0.06 

0.33 

0.23 

0 

+ n 

+ 


= 0.1539 < Co 


H 2 > max 

l< 2 <n 


1 - 




k=l,k^i 


-1 n 


A 




“ C,- Q;,+ I l<i<n 

J = 1 ^ V 


= max {2.7102, 0.8431, 0.6449} = 


= 0 . 2 , 
= 0.28, 
= 0.4, 

2.7102. 
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Figure 1: T = M. Numerical solution Xi{t) of system (4.1) for = 

(1.3,1.3,1.5). 

Let A 2 = 2.72, q = 0.7215, min{a)^} = 0.6, we have 

= min{1.2118,8.0234,1.3407} 

= 1.2118 > Til > ^ 0-7215 ^ ^ 2025. 

tty 0.6 

If —Cj G TZ~^, that is, 1 — Ci(t)fi{t) > 0,i = 1,2,3, then it is easy to verify that all 
conditions of Theorem 14.21 are satisfied. Therefore, the system in Example 4.1 has a unique 
positive almost periodic solution in the region B* = {(p| (p G B, Ai < (pi(t) < 2.72, t G T, i = 
1 , 2 ,..., n}, which is exponentially stable. 

Especially, if we take T = M or T = Z, then 1 — Ci{t)fi{t) > 0,i = 1,2,3. Hence, in this 
case, the continuous-time Nicholson’s blowflies model fll.2p and its discrete-time analogue 
have the same dynamical behaviors (see Figures 1-8). 


mm 

l< 2 <n 


1 - 


n 

E 

k=l,k^i 


ik 

ct 


20 














1.3 


Figure 2: T = 
(1.3,1.3,1.5). 


Figure 3: T = 
(1.3,1.3,1.5). 


1.2 - 
1.1 r 
1 1 

0.9 y 

><^ 1 
0.8 i 

1 

0.7 -i 

I 

I 

0.6 -1 
0.5 - 

0.4 I-^^^^^^^^^^- 

0 50 100 150 200 250 300 350 400 

time t 


M. Numerical solution X 2 it) of system (4.1) for {ipi{t),ip 2 it),ip 3 {t)) 



M. Numerical solution x^it) of system (4.1) for {ipi{t),(p 2 it),ip 3 {t)) 
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1.5 


Figure 5: T 
(1.2,1.2,2.3). 


Figure 6: T 
(1.2,1.2, 2.3). 



Figure 4: Continuous situation (T = M) : xi{t), X 2 (t), x^it). 



Z. Numerical solution xi(n) of system (4.1) for {ipi{n), ip 2 {n), = 



50 100 150 200 250 300 

time n 


Z. Numerical solution X 2 {n) of system (4.1) for {ipi{n), ip 2 {n), ip^^n)) = 
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Figure 7: T 
(1.2,1.2, 2.3). 



50 


100 


150 200 250 300 

time n 


Z. Numerical solution X 2 ,{n) of system (4.1) for ((^i(t), <y92(t),‘^ 3 (^)) = 



Figure 8: Discrete situation (T = R) : xi{n), X 2 {n), xz{n). 
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6 Conclusion 


In this paper, we proposed a new concept of almost periodic time scales, two new definitions 
of almost periodic functions on time scales and investigated some basic properties of them, 
which can unify the continuous and the discrete cases effectively. As an application, we 
obtain some sufficient conditions for the existence and exponential stability of positive almost 
periodic solutions for a class of Nicholson’s blowflies models on time scales. Our methods and 
results of this paper may be used to study almost periodicity of general dynamic equations 
on time scales. Besides, based on our this new concept of almost periodic time scales, one can 
further study the problems of pseudo almost periodic functions, pseudo almost automorphic 
functions and pseudo almost periodic set-valued functions on times as well as the problems 
of pseudo almost periodic, pseudo almost automorphic and pseudo almost periodic set-valued 
dynamic systems on times and so on. 
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